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In the slave boson mean field treatment of the t-J model, the ground state for small doping is 
a d-wave superconductor. A conventional superconductor has collective modes associated with the 
amplitude and phase of the pairing order parameter. Here, the hopping matrix element Xsj is also 
a mean field order parameter. We therefore expect that new collective modes will be introduced 
with these additional complex degrees of freedom. We compute the new collective modes and their 
spectral functions by numerically diagonalizing the matrix which describes the fluctuations about 
the mean field solution. We also show that the SU(2) gauge theory formulation allows us to predict 
and classify these collective modes. Indeed, the SU(2) formulation is essential in order to avoid 
spurious collective modes as the doping x goes to zero. The most important new collective modes 
are the 9 mode and the longitudinal and transverse <f> gauge modes. The 9 mode corresponds 
to fluctuations of the staggered flux phase and creates orbital current fluctuations. The 4> gauge 
modes correspond to out-of-phase fluctuations of the amplitudes of the pairing and the hopping 
matrix elements. We compute the neutron scattering cross-section which couples to the 9 mode and 
inelastic X-ray scattering cross-section which couples to the fluctuation in the real part of Xij- ln 
addition, we show that the latter fluctuation at (tt, n) may be coupled to the buckling phonon mode 
in the LTO phase of LSCO and may be detected optically. Experimental searches of these collective 
modes will serve as important tests of this line of attack on the high T c problem. 
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I. INTRODUCTION 

It is now a widely accepted view that the problem of high T c superconductivity is that of doping into a Mott 
insulator. The simplest model which captures the physics of the strong correlation inherent in this problem is the t-J 
model. The no-double-occupant constraint leads naturally to a gauge theory. The mean field decoupling of this theory 
is the formal language which describes Anderson's physical idea of a resonating valence bond (RVB). 1 The mean field 
theory indeed enjoys a number of successes. Notably it predicted the appearance of d-wave superconductivity and the 
existence of a spin gap phase in the underdoped region. 2-4 The properties of this spin gap phase are in remarkable 
agreement with those of the pseudogap phenomenology. Most striking among these is that the quasiparticle spectral 
function becomes a coherent peak with small spectral weight. This property emerges naturally out of the mean field 
picture of the condensation of bosons. 

Despite these successes, the theory has not enjoyed wide following because the gauge fluctuations are strong and the 
theory does not have a well controlled small expansion parameter, except formal ones such as the large N expansion. 
For example, Nayak has raised a number of objections on general grounds. However, we believe these general 
arguments have been adequately answered in the ensuing comments and discussions. 6-9 Here we briefly summarize 
our point of view. 

Initially, the gauge field has infinite coupling in order to enforce the constraint on a lattice scale. Upon integrating 
out high energy fermionic degrees of freedom, the coupling constant becomes of order unity. Then there is a chance that 
the mean field theory and fluctuations about it may be qualitatively correct at intermediate temperature scales. At low 
energy scales, nonperturbative effects related to compact gauge theories may come in, giving rise to the phenomenon 
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of confinement. This phenomenon is well known at half- filling. The mean field solution gives a 7r-flux state with 
Dirac spectrum centered about (±^,±1^). Coupling with compact U(l) gauge field leads to confinement and chiral 
symmetry breaking in particle physics language, which is equivalent to Neel ordering. 10 The idea is that with doping, 
the appearance of dissipation suppresses confinement, 11 and the Neel state is rapidly destroyed. Then the mean 
field description may be qualitatively correct beyond a small, but finite, doping concentration. At low temperatures 
characterized by the boson condensation scale xt, superconductivity emerges. In this theory the superconducting 
state is described as the Higgs phase associated with bose condensation. (We prefer to refer to this as the Higgs phase 
rather than bose condensation because the bose field is not gauge invariant.) However, here lies one of the significant 
failures of the weak coupling gauge theory description. As long as the gauge fluctuation is treated as Gaussian, the 
Ioffe-Larkin law holds and one predicts that the supcrfluid density p s (T) behaves as p s (T) fa ax — bx 2 T. The ax term 
agrees with experiment while the —bx 2 T term does not. 12,13 This failure is traced to the fact that in the Gaussian 
approximation, the current carried by the quasiparticlcs in the superconducting state is proportional to xvp. We 
believe this failure is a sign that nonperturbative effects again become important and confinement takes place, so that 
the low energy quasiparticles near the nodes behave like BCS quasiparticlcs which carry the full current vp- For the 
pure gauge theory, the confinement always occurs in (2+l)D however small the coupling constant of the gauge field 
is. When the bose field is coupled to the gauge field, i.e., Higgs model, it is known that the confinement phase are 
smoothly connected to each other. 14 Recently, we showed that in (2+l)D Higgs model the phase should be considered 
confined everywhere. 15 Thus the appearance of confinement at low energy scale is not surprising, but the dichotomy 
between the success of the bose condensation phenomenology on the one hand (in explaining the quasiparticle spectral 
weight mentioned earlier) and the confinement physics needed to give the proper quasiparticle currents on the other, 
is in our opinion one of the most profound questions facing the field today. 

To summarize, the current status of the gauge theory approach to the t-J model is the following. Within mean field 
theory, a certain saddle-point solution may be regarded as the "mother state" at some intermediate energy scale. The 
most promising candidate for this "mother state" appears to be the staggered flux state, which is the saddle point 
for the SU(2) formulation of the t-J model. This saddle point is an unstable fixed point due to gauge fluctuations. It 
flows to the Neel state at or near half filling and to the superconducting ground state at finite doping. This picture 
is in accordance with a recent insightful comment by Anderson, 16 except that we provide a concrete description of 
the unstable fixed point. As discussed above, the issues raised by this picture are profound and probably difficult to 
treat analytically. In the past few years, we have focused on trying to substantiate this point of view by numerical 
methods and by predicting new experiments. On the numerical front, the study of projected wavefunctions have 
yielded remarkable insights. Ivanov et al. 17 reported a finding which is surprising except from the SU(2) gauge theory 
point of view. They discovered that upon Gutzwillcr projection of a BCS c?-wave state, the current-current correlation 
shows a staggered pattern, as expected for the staggered flux state. Thus there is strong numerical evidence that 
the staggered flux state and the d-w&ve superconducting state are intimately related. On the experimental front, 
predictions have been made that the staggered flux state may be stabilized in the vortex core and ways to measure 
the effect of the orbital currents have been proposed. 18 In this paper we continue work in this direction and ask the 
question: is the superconducting ground state that emerges out of the gauge theory completely conventional, or are 
there detectable consequences of its unconventional origin? 

The short answer to this question is that gauge theory predicts new collective modes in the superconducting ground 
state which have experimental consequences. Unlike conventional BCS theory, where the only order parameter is the 
complex pairing order parameter A, here the hopping matrix element % also functions as an order parameter. Thus 
it is natural to expect new collective modes. This problem is formulated generally in Sections II and III. In Section 
IV we show how the gauge theory allows us to predict the low lying collective modes. The new modes are the 6 mode 
and the transverse and longitudinal cj> modes. The 6 mode generates staggered orbital current fluctuations and is 
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related to the current fluctuations found in the projected wavefunctions. However, while the projected wavefunctions 
give only equal time correlation, here we obtain the full dynamical spectral function. The <f> gauge modes are new 
excitations related to amplitude fluctuations of \ an d A. In Sections V and VI we present numerical results of the 
spectral functions and describe experiments which may couple to them. 

The approximation employed in this paper is to treat the system in the Higgs phase by adopting the radial gauge for 
the bosons. Gauge fluctuations are treated at the Gaussian level and nonanalytic corrections, such as instantons which 
may lead to confinement, are neglected. This can be justfied in terms of the l/7V-expansion with N being the number 
of fermion species. In the present case, N = 2 and one criterion for the validity of this Gaussian approximation 
is that the magnitude of the order parameter fluctuation is less than that of the mean field value. As is evident 
from the discussion in II and III, the magnitude of the fluctuation diverges as x — > because the action for the 
SU(2) local gauge transformation vanishes in this limit. Hence the Gaussian approximation breaks down near x = 0, 
where the confinement physics is of vital importance because the nonperturbative configurations such as instantons 
contribute appreciably. Below we will show mainly the numerical results for x = 0.1. We have estimted numerically 
the magnitude of the fluctuation of the order parameters, e.g., < (Sxx) 2 > at zero temperature, and found that it is 
of the order of unity at x = 0.1. Therefore the gaussian approximation is marginaly justified in this case. As another 
test, recently Honerkamp and Lee 19 have computed the shift in A using a free energy which includes the gaussian 
fluctuaions of the 9 mode. They found a reduction to 55 % of the mean field saddle point value for x=.06 but for 
x=.l the reduction is only 80 %. Thus the effect of the gaussian fluctuation is relatively under control for increasing 
doping. Furthermore, since the collective modes appear at fairly high energy (~ J) for experimentally relevant doping 
(x s=s 0.1), it is possible that this approximation is valid at this energy scale, but the ultimate test must come from 
experiment. 



II. U(l) AND SU(2) FORMULATIONS OF THE T-J MODEL 



We begin with the familiar U(l) formulation of the t-J model 20 

H = j(Si-Sj- jn.nA - h jC \ a c ]a . (1) 

<i,j> ^ 

The constraint of no double occupation is enforced by writing 

4 = fl^ (2) 

and imposing the condition J2<j fLf™ + ^j&i = 1> which m turn is enforced with a Lagrangian multiplier Aj. The 
Heisenberg exchange term is written in terms of fi a 21 

Si ■ Sj = —-^fivfjafjpfiP 

\ (44 - 44) C/ii/it - fnh) (3) 



+ 

We write 



mrij = (i-b\b i ){i-b]b j ). (4) 

Then Si ■ Sj — jriiUj can be written in terms of the first two terms of Eq. (3) plus quadratic terms, provided we ignore 
the nearest-neighbor hole-hole interaction jb\bibjbj. We then decouple the exchange term in both the particle-hole 
and particle-particle channels via the Hubbard-Stratonovich transformation. By introducing the SU(2) doublets 22 ' 23 
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J a J V /»t 

the partition function is written in the compact form 



(6) 



Z = y DQD&DbDXDU cxp ^- drL^j 



where 



Li = ^ E Tr [4^] + i E (C^^ + cc 

<ij> <ij>,a 
+ E^r-^)./W 

+ ^6*(<9 r — iAi + /u)6i 

-J^tyb^flfj,, (7) 



with representing fermion hopping and Ay representing fermion pairing. In Eq. (7) Jy = J/4 but in the literature 
J has sometimes been taken to be 3 J/8. The latter has the advantage that the mean field equation reproduces that 
which is obtained by the Feynman variational principle, 24 but these differences are well within the uncertainties of 
the mean field theory ( In the mean field theory \ij = Y. a < fLfjv > and A y =< Mfjl ~ /a/jT >•) 

Affleck et al. 22 pointed out that the t-J model at half-filling obeys an exact SU(2) symmetry in the functional integral 
formulation. The SU(2) doublet in Eq. (5) expresses the physical idea that a physical up-spin can be represented 
by the presence of an up-spin fermion, or the absence of a down-spin fermion, once the constraint is imposed. Wen 
and Lee 25 proposed a formulation which obeys the SU(2) symmetry even away from half-filling. The SU(2) and the 
original U(l) formulation arc equally exact, but once approximations are introduced, the SU(2) formulation has the 
advantage that the zero doping limit x — > can be smoothly taken. We shall see an example of this in the collective 
mode spectrum described below. In the SU(2) formulation a doublet of bosons is introduced 

fc-(fc). (9) 

The physical Hilbert space is the SU(2) singlet subspace. The electron operator is an SU(2) singlet formed out of the 
fermion and boson doublets 

Cia = 75 h ^ ia (10) 

and three Lagrangian multipliers a^, I = 1,2,3 are needed to project to the SU(2) singlet subspace and impose the 
constraints 



o*t^» + htrhi = . (11) 



Now the partition function Z is given by 



Z = J D&D&DhDakDcftDcftDUexp ^- drL^j 
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with the Lagrangian taking the form 



J 



J 
2 



(12) 



As pointed out in Ref. (24), Eq. (12) is closely related to the U(l) Lagrangian Eq. (7) if we transform to the radial 
gauge, i.e, we write 



hi = Qi 



where bi is complex and gi is an SU(2) matrix parametrized by 



Zi2 Z*\ 



where 



and 



_ -±i 1 

Zn = e M, e 1 2 cos ■ 



Z12 = e Ml e 4 2 sin— . 



(13) 



(14) 



(15) 



(16) 



The angle on in zn and Zi2 is the overall phase which is redundant and can be absorbed in the phase of bi. 
An important feature of Eq. (12) is that L 2 is invariant under the SU(2) gauge transformation 

hi = g\hi 
Uij = Uj 

and 

«« t£ = glaii^gi + ig\ {d T gi) ■ 



(17) 
(18) 
(19) 

(20) 



Starting from Eq. (12) and making the above gauge transformation, the partition function is integrated over bi and 
Qi instead of hi and the Lagrangian takes the form 

J 



<ij> 



<ij>,cr 



+ lY,$l(dr-ia i oiT*)$ i 

i,<7 

+ E 6 * & - ia oi + ») b * 

i 



(21) 
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We have removed the tilde from Uij, <&i a , fi a , Oq because these are integration variables. Note that gi has disappeard 
from the actions and Eq. (21) is the same as the U(l) Lagrangian L\, with the exception that Uj is now replaced 
by Uj — Uj/2, Xi becomes and, most importantly, two additional integrals and aj^ coupling to the fermions 
appear. We note that in the limit of zero doping, thanks to these additional gauge fields, Eq. (21) manifestly invariant 
under SU(2) transformation, whereas Eq.(7) is not. This will have important consequences when we consider mean 
field approximation and small fluctuations, in that Eq. (21) will have smooth x — > limit while Eq. (7) does not. 



III. MEAN FIELD THEORY AND COLLECTIVE MODES 



We now consider the mean field treatment of Eq. (21) and the quadratic fluctuations about the mean field, which 
will yield the collective modes. We work in the radial gauge, where 6, is considered real without loss of generality. As 
discussed after Eq. (16), the phase of hi and are redundant and one of them can be chosen as zero. We will discuss 
the alternative choice later, but here we choose hi to be real. The saddle point solutions are 

bi = r 
Xij = Xo 

A ii<+ „ = Aor/p or A itj = Aoi-l) 1 ^ 

*(& = (0,0,Ao) (22) 

where /i = x or y and r\ x — l,r) y — — 1 correspond to (i-wave pairing of the fermions. The saddle point corresponds to 
a physical (i-wave superconductor, as the order parameter (cjfCj+^j — c^Cj+^-f) = r^lS.^^ is nonzero. The mean field 
fermionic action is 



E { A 



where \k = ^Jxo(cosk x + cosfc y ), A^ = 2JAo(cosfc x — cosk y ), = 2trQ(cosk x + cosfc y ) and ^ = — Xk — tk — Ao- 
We write the small expansion about these saddle points as 

bi = r (l + SRi) 
= Ufj + 5U l3 

<4 = ($ a L tali, -* x o + Sa^) (24) 
where the mean field U d describes d-wave pairing, 

U? ii+I1 - -Xor 3 + Aq^t 1 . (25) 



The fluctuation is expanded as 



where r° = /. Note SU X and SUy,a = 1 to 3 are real variables, while for a = are purely imaginary, and together 
make up a total of 8 degrees of freedom. These correspond to complex hopping Xi,i+n an d pairing A iji+A1 in the x 
and y directions. 

By setting the linear terms in the small expansion of the free energy to zero, we obtain the standard saddle point 
equations A.l to A. 3. The second order deviation is described by a 12 x 12 matrix, where the variables are 6U£, SR, 
and 5qq. The details are given in the appendix. As it stands the matrix is not hcrmitian. On the other hand, if we 
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consider ia as variables, the matrix is hermitian for w n = but has negative eigenvalues, i.e., it corresponds to the 
saddle point of the free energy in the unprojected Hilbert space. In order to obtain positive eigenvalues, it is necessary 
to first integrate out Sa^ in order to project to the physical subspace. Since this is a Gaussian integration, we may 
equivalently consider (5U®, 5R) as the physical degrees of freedom, and solve for the local Sa^ for each configuration. 
This is the idea behind the cr-model approach in rcf. (24) where large fluctuations in SUij are considered. The present 
work should be considered the low-temperature limit of the cr-modcl. 

At this point we proceed numerically, evaluate the 12 x 12 matrix and integrate out the <5oq fields. The remaining 
quadratic form gives a 9 x 9 matrix with 9 eigenvalues. As expected, there is a soft mode associated with the phase of 
the pairing order parameter. In addition, we find a number of soft modes in the small x limit. Before presenting the 
numerical results, we show how the SU(2) symmetry allows us to predict and classify all the soft degrees of freedom. 



IV. SU(2) CLASSIFICATION OF SOFT MODES 



In this section we make use of the SU(2) gauge symmetry to classify the soft modes. The basic idea is the following. 
As the temperature is decreased, the SU(2) symmetry is broken via a series of symmetry breaking at the mean field 
level. For small x, SU(2) is first broken down to U(l) at a temperature scale of order J/2 to the staggered flux (s-flux) 
state. 25 At a lower temperature of order xt, the bosons condense and the gauge symmetry is broken completely. Of 
course, a local gauge symmetry cannot be broken, but the mean field description is still a useful starting point to 
describe the low-lying collective excitations, which are physical. A familiar example is the pairing order parameter of 
superconductivity, which breaks the local U(l) gauge symmetry associated with the E& M gauge field at the mean 
field level. While strictly speaking, this order parameter is not gauge invariant, it is a useful starting point which 
leads to the correct description of the gauge field via the Anderson-Higgs mechanism. 

It is useful to distinguish between two kinds of symmetry breaking as x — > 0. First, at x — the mean field solution 
is the 7r-flux state, i.e., xo = Ao- This state has full SU(2) symmetry, which is broken down to U(l) in the s-fiux state, 
where xo Ao. We shall refer to the remaining symmetry as the residual U(l) symmetry. As x becomes nonzero, 
Ao/xo deviates from unity rather rapidly and we shall focus our attention on the zero modes due to the residual U(l) 
symmetry. Secondly, the bosonic degrees of freedom appear at x ^ 0. Boson condensation breaks the residual U(l) 
completely below an energy scale of order xt. The zero modes then acquire a finite energy gap which is the subject 
of our analysis. 

Starting from the d-w&ve superconductor mean field solution described in the last section, we expect the soft modes 
to involve small fluctuations of the boson about the radial gauge which can be parametrized by the SU(2) matrix gi 
such that hi = gi (^^j . In addition, we include phase fluctuation of the Uij matrix which wc paramaterize by 

Un = Uf j <i ia »- r (27) 

where afj, £ = 1, 2, 3 are three gauge fields living on spatial links. Since SU(2) has been broken down to U(l), only 
one out of three gauge fields remain soft in the s-flux state. 27 ' 28 To visualize this, it is convenient to make a gauge 
transformation using Eqs. (17, 18, 19) to the s-flux order parameter 

U?f = w\uf w 3 (28) 

where 



Wj — exp 

and 



i(_l)j.+j»| T ij (29) 
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= -Xor 3 - iA (-iy*+iy 

= -At 3 cxp(z(-1)^ +j «$ t 3 ) (30) 

where %o = Acos$o and A = ^4sin<J> - Equation (30) represents fermion hopping with a complex matrix element 
such that a flux 4$ threads the lattice plaquettes in a staggered manner. At the same time the boson is transformed 
to 

-3 (-«-?>-.) (31) 

and the mean field a becomes a,Q F = (0, iA (— l) Jx+Jy , 0). We note that Ufj F describes a semiconductor band with 
nodes at (7r/2, 7r/2). If hi were to remain as hi = r ( J and ajj ^ , we would have described an s-flux state 



with small fermion pockets. Instead, hi and ag F are rotated such that Oq F couples to pair fields fj^fj^ and the 
resulting state is gauge equivalent to the d-w&ve superconductor that we started out with. The advantage of the 
s-flux representation is that Ufj F is proportional to T3 and is invariant under T3 rotation. Thus the residual U(l) 
symmetry is apparent. We expect the soft modes to be described by 

W = Uf F e^ (32) 

and 

0) ( 33 ) 

where h\ is close to hf F , i.e., g\ is parametrized by Eq. (14) with 9 close to | and <f) close to (— l) lx+% y^. In Eq. (32) 
we have ignored the ajj and afj gauge fields as they have been pushed to finite frequencies by the Anderson-Higgs 
mechanism. 

In order to visualize the different gauge choices, it is useful to introduce the local quantization axis 

Jj = zjrzi = (sm9iCOS(j)i,sm6ism(fii,cos6i) . (34) 

Note that Jj is independent of the overall phase a.{. In the s-flux representation the quantization axis has been 
rotated with Wi given by Eq.(29) to point along the ±y-axis in a staggered fashion. Small fluctuations correspond to 
59 deviation from the equator and S(f> in the azimuthal angle. This is illustrated in Fig.l. 



8 



staggered flux Q Q 




staggered flux Q Q 



FIG. 1. The quantization axis I in the SU(2) gauge theory. The north and south poles correspond to the staggered flux phases 
with shifted orbital current patterns. All points on the equators are equivalent and correspond to the d-wave superconductor. 
In the superconducting state one particular direction is chosen on the equator. There are two important collective modes. The 
modes correspond to fluctuations in the polar angle 89 and the <f> gauge mode to a spatially varying fluctuation in 8<j>. 

It is useful to rotate the configuration specified by Eqs.(32) and (33) back to the radial gauge. We obtain 

Ut^tfufife^j. . (35) 

The advantage of the radial gauge is two-fold. The electron operator Ci a — rofi a and we can consider Ufj as the 
effective Hamiltonian for electron quasiparticles. At the same time, we can now make contact with the fluctuation of 
the Uij matrix in the last section and interpret the numerical results. 

Equation (35) can be explicitly evaluated for arbitrary 0, (f> and a and a?-, resulting in 



U' 



(36) 



The overall phase on enters the effective hopping U22 and effective pairing Uf 2 in the expected way, and Ufj is the 
a>i = limit given by ref. [27]. 



~Xij 



T 3 COS • 



i(-iy 



+ (-1) 



J*+% T 2 gin . 



v COS • 



'} 1 • 

— — t sin ■ 



(37) 



where 



Xij = A cos $ij 
Aij = Arjj-ism&ij 



(38) 
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and 



va = 



— ^-a 3 - 
2 °« 



(39) 



Note that the only dependence on 0;, <f>j is via the gauge invariant combination v^, which has the interpretation of 
the gauge current. Furthermore, for 9i = 9j = §, we see from Eqs.(36) and (37) that Xij^ l<<ai ~ a ^ and Ajje~'( ai+a ^ 
play the role of the effective hopping and pairing parameters. Thus fluctuations in <j>i leading to nonzero Vij means a 
fluctuation in the amplitude of Xij and A i3 in such a way that \x,ij\ 2 + \^ij\ 2 = A 2 is fixed. This is shown in Fig. 2. 




X 



v 



FIG. 2. Geometrical interpretation of the fluctuation of Vij. The angle Xij [given by Eq.(38)] is modulated around the flux 
<I>o in a staggered manner, in such a way that the hopping amplitude 4^ and the pairing amplitude A^- are modulated, also 
in a staggered manner. 



To make contact with the collective modes, we write 9i = \ + S9i and expand Eqs. (36,37) to first order in 59i, Vi 



and on. 



Ufi - U* = -xor 2 i l)l * +lV {59, - 69j) + iA oVft { W + S9 3 ) 



+ (xot 1 + A ^r 3 ) 



+ ^Ao^r 2 (ai + ctj) 



3 . 



(40) 



where j = i + jj,, jj, = x, y. 

Equation (40) is a main result of this section, as it allows us to interpret the collective modes when compared with 
Eq.(26). It predicts the location in momentum space of the soft collective modes and gives the eigenvectors. We 
classify the modes as follows: 

uniform a : on w a v , 



staggered a s (-l) lx+ly 

uniform 9 : 9i rts 9 U , 
staggered 9 : 9i w S (-1 )'*+'" 



(41) 
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We consider a u , a s , 9 U , 9 S and Uj^+^Qu — x, y) as six slowly varying variables and the first four variables are defined 
as follows 



(aj + otj) = a u 



-{ai-aj) = (-l)^ +i "a s 
- (0j + 8j) = 9 U 

\Wi-9i) =(-iy^e s (42) 

Substitution into Eq.(40) and comparison with Eq.(26) show that at q = 0, we can identify the following collective 
modes and their eigenvectors. 



Goldstone mode: a„ , XJl — Uy 



(43) 



internal phase mode: 9 S , U% + Uy 
At q — (tt, tt) we pick out the coefficients of (— 1)*»>+ Ii » and identify the following modes and corresponding eigenvectors. 



9 mode: 9 U , U^-U y 



°-U° 

x y 

<j) gauge mode: v iti+fi , xoUl + A Q L^ , xoUy - A U* (44) 
X mode: a s , U° + Uy 

The notation aU" + bU@ means that each component of the eigenvector is (6U£ = a,SU^ = b), etc. in Eq.(26) The 
nature of the collective modes are readily identified from their eigenvectors. Two of these modes were known before. 
The Goldstone mode a u is the standard one associated with the phase a u of the superconducting order parameter 
A M e~ l2a ". The 9 S mode corresponds to the out-of-phase oscillation of the phase of the superconducting order pa- 
rameter in the x and y directions, lA^le -12 ^ 31 and |Aj / |e _l2Vy , such that 9 S — \{(p x — <py)- This is a property of any 
d-wave superconductor and is labelled the internal phase mode. 

The a s mode corresponds to the fluctuation in the phase of Xij- hi BCS theory the hopping term is fixed and not 
allowed to fluctuate. So this is a new degree of freedom special to the gauge theory. The slowly varying phase of \ij 
plays the role of the spatial component of the gauge field in the U(l) gauge theory. 20 This appears as a collective 
mode at q = (tt, tt). We shall see that just as in the U(l) theory, the a s mode plays the crucial role in producing the 
correct answer for the superfluid stiffness p s = x. 

The new modes that are of greatest interest to us in this paper are the 9 mode and the <j> gauge modes. These will 
be discussed in greater detail later. 



staggered 8-mode 



F, 0.3- 



staggered <j>-mode 



// 
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FIG. 3. The eigenvalues at q = (0, 0) and u> = of the 9x9 matrix which describes fluctuations about the mean field as a 
function of doping x. Shown are the eigenvalues which vanish as x — > 0. The eigenvalues are given in units of J. In addition 
to the Goldstone mode corresponding to superconducting phase fluctuations which has zero eigenvalue for all x (not shown), 
we find the staggered 8 mode [called internal phase mode in Eq.(43)] and a continuation of the transverse 0-gauge mode from 
q = (7r,7r). 

Now we compare our analysis with numerical results described in the last section. In Figs. 3 and 4 we plot all the 
eigenvalues of the 9x9 matrix which vanishes at x — as a function of x at uj = and at q = and q = (tt,tt), 
respectively. They can all be identified with our classification. At q = the Goldstone mode has zero eigenvalue 
for all x as expected. The internal phase mode rises rapidly with increasing x. In addition, we find a mode with an 
eigenvector that corresponds to the continuation of the transverse <p gauge mode we find at q — (7r,7r) to q = 0. 



staggered a-mode 



> gauge modes 



uniform 9-mode 



FIG. 4. Same as Fig. 3 but at q = (7r,7r). We find the uniform 9 mode, two degenerate <j) modes, and a staggered a mode 
[called the \ mode in Eq.(44)]. 

At q = (n, 7r) we find the 9 mode and the two cj> gauge modes (transverse and longitudinal) which are degenerate. 
In addition we find the x mode. The vanishing of the eigenvalues of all these modes at x = is a consequence of 
the SU(2) symmetry. However, as we discussed before, even if |x| ^ |A|, as long as x = the 9 mode still has zero 
eigenvalue and the <fr gauge modes are protected by the residual U(l) symmetry. These modes are gapped only by 
the boson condensation. 

We should mention that in addition to these soft modes, we found two unstable eigenstates for small x. If we 
use in Eq. (21) t/J = 16/3, i.e., t/J = 2 and t' = — 0.5t, which generates a realistic Fermi surface, the eigenvalue 
corresponding to r fluctuation is negative for x < 0.084. Since Tq corresponds to hole density, this signals an instability 
to phase separation. This instability is easily suppressed by long range interaction and is largely decoupled from the 
modes of interest. A second instability occurs for x < x c near q = (0, w) and (n, 0). The eigenvector corresponds to 
amplitude fluctuations in \\\ 2 + |A| 2 which we identify as instability to columnar dimer formation. This instability 
is well known for x = and can be suppressed by bi-quadratic terms. 28 For the parameters mentioned above, we 
find x c — 0.087. Most of our detailed numerical results are for x > x c . Furthermore these two modes are mostly 
decoupled from the low lying modes in Eq.(41), and does not disturb their behavior in the limit x — > discussed in 
the followings. 

We next discuss in greater detail the 9 and the <j> gauge modes. As seen in Fig.l, the 9 mode is a fluctuation of 
i" towards the north and south poles, and describes the admixture of the s-flux phase. From Eq.(37) we see that 6 U 
generates a staggered imaginary part to the hopping matrix elements which will produce staggered orbital currents. 
Thus the physical manifestation of the 6 mode is staggered orbital current fluctuations. The softness of the 9 mode 
means that it is readily excited by thermal or quantum fluctuations. This is clearly related to the strong staggered 
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orbital current fluctuations found in the projected d-wave superconductor wavcfunction. 17 The energy cost is low 
because the s-flux and the d-wave superconductor are almost degenerate in energy. The energy difference arises only 
because clq is different in the two states. At the mean field level, the energy difference comes from the Fermi pockets 
and is proportional to x 2 . However, we find that after integrating out a , the energy cost increases, apparently due 
to the enforcement of the constraint and the eigenvalue shown in Fig. 3 is linear in x. We have also computed the 
eigenvalues for finite u>. It should be noted that the 9x9 matrix is not Hermitian for finite w m (Matsubara frequency). 
Then we make the analytic continuation as ioj m — > u + id with infinitestimal S > 0. This S can be neglected outside 
of the particle-hole continuum, and in this case the 9x9 matrix is Hermitian and its eigenvalues are real. We will 
present below the results for (q, u) outside of the particle-hole continuum. We found that the eigenvalue can be fitted 
by ax — |o> 2 — fk 2 for small uj and k = q — (tt, tt). The inverse x dependence of the coefficient might be surprising 
but the to 2 and k region shrinks as x becomes small, so that the x = limit is smoothly attained. Furthermore, the 
negative coefficient of k 2 indicates that the eigenvalue has a local maximum at q — (tt,tt). A plausible dependence is 
y 7 (ax) 2 — b'uj 2 — c'k 2 , which is reduced to the above form in the limit \u>\, |k| << x. Back to the Matsubara frequency, 
the suggested effective action for the 9 mode is 

Se=J2 \/W 2 +H-^|«ft^)| 2 (45) 

k,ui m 

in the region \ui m \, |k| << x. The small energy gap ax leads to a strong spectral weight. This is confirmed by direct 
numerical calculations in the next section. Here it is interesting to compare these results with the SU(2) formalism 
with those in the U(l) where only a 3 , is integrated over. The suggested action in this case is 

So {1) = E + b "" 2 ™ + C " xk2 } m,Urn)\ 2 . (46) 

Therefore it is evident that the SU(2) symmetry leads to quite different x-dependence in the limit x — > 0. 

To investigate the <j) gauge modes we first discuss x = 0. The mean field solution is the 7r-flux phase where the 
fermions obey the Dirac spectra with nodes at (±|,±|). After integrating out the fermions, it is known that the 
effective gauge field action is purely transverse and given by 30 

(k k \ 
S ' iV kT~ ) M fc > u m )a„(k, -u) m ) (47) 

where a M is the continuum version of a^i+^j and k is measured relative to (7r, tt). We confirm this by computing 



the eigenvalues for x = at finite q = (1.037r,7r). The transverse mode behaves as \/k 2 + uj^ as expected while 
the longitudinal mode is exactly zero for all ui m . This is because the longitudinal mode is pure gauge (a M = V^) 
and is not a real degree of freedom. In constrast, if we worked with the U(l) formulation [Eq.(7)], we find that the 
longitudinal mode has eigenvalue close to Thus the addition of the aj and a 2 , in the SU(2) formulation is crucial. 
Otherwise we would have gotten a spurious collective mode. This is a dramatic illustration of the advantage of the 
SU(2) formulation, if one is interested in obtaining meaningful results which are smoothly connected to the undoped 
case. For finite x we see from Fig. 3 that the eigenvalue increases linearly with x for small x. A reasonable approximate 
for the transverse mode is [for k = q — (tt, tt)] 



^transverse 

where 



= E {^Q\J ( a v x ) 2 + fe2 + w m ) (<W ~ V n( k ^m)v v (k,u; m ) (48) 



>n = 2 V ^ - °m (49) 
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is the continuum limit of 1^,.;+^ and ceo parameterizes the spectral weight. Thus we expect the transverse mode to 
show an edge singularity starting at uo = (a^x) 2 + k 2 . In Fig. 5 we show the eigenvalues of the </> gauge modes as 
function of (q x — ir)/it at q y = 7r and u) m = 0, which shows the expected dependence on q x — ir from eq.(48) for 
transverse mode. 
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FIG. 5. The eigenvalue (in units of J) of the transverse and longitudinal c/> gauge mode as a function of k x /ir — (q x — 
with q y — 7r and u> m = for x = 0.004. They show the expected behavior for small x and q near (it, it) because when x — 
and ui m = 0, the transverse mode should be linear in \k\ and the longitudinal mode should be zero for all q and ui m . 



At finite x the longitudinal mode becomes a real degree of freedom due to the breaking of the residual U(l) symmetry. 
At w m = and q = (n, n) the eigenvalue is a^x, degenerate with the transverse mode. For finite k x = q x — 7r, the 
eigenvalue shows a slow decrease and then saturates. This is shown in Fig. 5. This is to be expected as the eigenvalue 
is zero for all u) m and q for x = 0. As for the spectral function, we just invert the 12 x 12 matrix numerically with 
uj m —> uj + i5, and we do not have to worry about the non-Hcrmitian nature of the matrix, which will be shown in 
the next section. It is worth noting that the longitudinal gauge mode is the only mode with significant coupling to 
5R and therefore to density fluctuations. 

We next discuss the Goldstone mode. We expect the eigenvalue to be of the form p s (Va) 2 , where a is the continuum 
limit of a slowly varying a u . We indeed verify that upon diagonalizing the 9x9 matrix described in the last section, 
a soft eigenvalue emerges which is linear in q 2 . Furthermore the coefficient p s is proportional to x. This is to be 
expected in our SU(2) formulation, as x = is an insulator. It is interesting to analyze how this result emerges. The 
point is that the coefficient of (Va u ) 2 in the free energy computed using Eqs. (21,40) is pf, which corresponds to the 
superfluid density of a nearly half- filled conventional superconductor, and pj ~ 1. How do we obtain p s = x ? In the 
numerics we find that the p s q 2 mode is a coupled mode between a u and a s - Recall that a s is related to the phase 
of Xij and may therefore be identified with the U(l) gauge field a in the U(l) formulation. 20 In that formulation the 
free energy in the superconducting (Bose condensed) state takes the form 



F = p f (Vc 



Pb\a\ 



(50) 



where pi, — x and pf is the fermion contribution to the superfluid density, which is of order unity. Upon minimizing 
F with respect to a, we arrive at the Ioffe-Larkin formula 30 F — p s (Va„) 2 where 



PbPf 
Pb + Pf 



(51) 



The screening of Va u by the gauge field a converts the fermion response pf to the physical response p s ~ x. When 
we express Eq.(50) in matrix form we see that the free energy has diagonal contributions p/(Va M ) 2 and (pf + pb)a 2 
and an off-diagonal term — pf Va„ • a. When we examine the 9x9 matrix, we find that for small q, a u is only coupled 
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to a s and the form of the 2x2 sub-matrix is just that given by the above discussion if we identify a with a s . Upon 
diagonalizing the 2x2 matrix, the soft mode with eigenvalue p s q 2 emerges. This simply confirms that our SU(2) 
formulation contains the same U(l) gauge field and the same screening mechanism is at work. One consequence of 
Ioffe-Larkin screening is that at finite temperature p s (T) = p s (0) — Cx 2 T, i.e., the coefficient of the linear T term is 
proportional to x 2 . We have also verified numerically that this is the case by computing p s at finite T. Experimentally 
there is strong evidence that the coefficient of the linear T term is independent of x. 31 The fact that p s is proportional 
to x can be seen more readily if we associate the phase Qj with bi in Eq.(13) instead of with as we have done so far. 
Then it is clear that for static a,, only the last term in Eq.(12) depends on a, — ctj and the free energy change must 
be proportional to 7q w x. Furthermore, the coupling to fermion excitations is proportional to x, so that the response 
to thermal excitations of quasiparticles is proportional to x 2 . Thus loffe-Larkin and our fluctuation theory are in 
disagreement with experiment. We believe this is an indication that the fermions and bosons are confined to become 
electrons in the superconducting state and that confinement physics is beyond the Gaussian fluctuation considered 
here. 

V. NUMERICAL RESULTS FOR THE SPECTRAL FUNCTIONS OF THE MODE AND THE GAUGE 

MODES 

As shown in Appendix A, the collective fluctuations about the saddle point is described by the quadratic form 

S eS = X a (M^ + M^)X l3 , a = 1...12 (52) 

where X a = (SU^, ...,5R, aj, ap, a§) are the 12 degrees of freedom and M£p(q, iuo m ) are made up of fermion bubbles 
computed in Appendix A. Here we study the spectral functions 

S A (q, uj) = ImC A (q, iw m — > w + iff) (53) 

where 

C A (q, iw m ) = (A(q, iu m )A(q, -iu m )) (54) 

and we focus on three cases, the mode (A = 9 U ), the transverse and longitudinal <j> gauge modes (A = Vis +fl ). Using 
Eq.(44) these are readily expressed in terms of 

{X a X p ) = (M J + M F )~ 1 (55) 

and computed numerically. We use parameters in Eq. (21) t/J = 16/3, i.e., i/J = 2 and t'/t = —0.5, which gives 
realistic fermion bandstructure and we show results for x = 0.1. The mean field parameters are xo — 0.376 and 
A = 0.255. Note that the maximum energy gap is E g = A fe= ( T Tr ) = 4A J We measure energy in units of J and the 
lattice constant is set to unity. 

It is important to note that we define the correlators in terms of the eigenvectors [(Eq.44)] which are the eigenvectors 
fo r q = (7r, 7r) and uj m — 0. Away from (tt, it), the overlap with the true eigenmode is modified and other modes may 
mix in. Furthermore, for the 0-gaugc modes, the assignment in Eq.(44) corresponds to "polarization" in the x and y 
directions, and are the appropriate longitudinal and transverse modes only for q along the (q x , tt) direction. Thus the 
numerical results shown here should be viewed as providing a guide for the behavior of the modes near (it, tt). In the 
next section we will compute correlation functions which are experimentally observable in the entire Brillouin zone. 
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X Axis Title 



FIG. 6. Spectral function for the 6 mode for q = (q x ,q y ) from to (n.n). A contour plot and a bird's eye view are shown. 
Frequency ui is in units of J(= | J). 

Figure 6 shows the results for the mode Se(q, lo) with q along the diagonal q x = q y in a contour plot and also in 
a bird's eye view. We see that the spectral function shows a strong peak near uj « 1.1J (J is defined after Eq.(8) and 
is suggested to take the value J = |J) which is strongly localized near q — (ir,n). The spectral function disperses 
rapidly upwards as q deviates from (it, it). The very large peak height can be anticipated from the approximate form 
given by Eq.(45) because the small value of the gap ax gives a large spectral weight upon inversion of the effective 
action. The strong and narrow peak is responsible for the orbital current fluctuations in the superconducting ground 
state. 17 




FIG. 7. Spectral function for the transverse cf> gauge mode for q — (q x ,n), i.e., from (0,7r) to (ir,iv). A contour plot and a 
bird's eye view are shown. Frequency u is in units of J(— | J). 

Figures 7 and 8 show the transverse and longitudinal <p gauge modes along q = (q x , it). The operators are chosen 
to be 

S(j) t = cos ^oSAy — sin &o5Xy (56) 
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5<f>t = cos ® 5A X + sin $ $Xx (57) 

where 

S X» = ^Y1 fiofi+Ha + c - c -) - XO (58) 
a 

SA H = \{f^fi+ixi ~ hlfi+rf + c-c) - Ao^ (59) 

are the fluctuations of the real part of \ij an d Ay, respectively. Since Ao are real, these correspond to amplitude 
fluctuations. These correlators are readily related to the correlations involving and [/*, respectively. Equations 
(56,57) are simply the resolution of the fluctuations illustrated in Fig. 2 into its components along the vertical and 
horizontal axes. Near (jr, ir) the two modes are degenerate and show a peak around uj — 1.5 J. The lincshapc is shown 
in Fig. 9. The mode is damped by particle- hole excitations and a sharp feature appears as the frequency drops below 
the particle-hole continuum. (Recall that the particle-hole continuum is set by the scale u> = 2E g w 2J with our 
parameters.) The transverse mode frequency remains low but loses spectral weight as q goes away from tt, it towards 
(0, 7r). On the other hand, the longitudinal mode disperses upwards and gains in strength. 




FIG. 8. Contour plot of the spectral function for the longitudinal cf> gauge mode for q = (q x ,iv), i.e., from (0,7r) to (ir,n). 
Frequency uj is in units of J(= | J). 
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FIG. 9. The lineshape for the c/> mode at q = (7T,7r). Frequency u is in units of J(= 
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Finally we discuss another important mode which is closely related to the <f> gauge mode. It is the fluctuation of 
the amplitude = \J\xi.i+^\ 2 + lA^+^j 2 and is parametrized by 



SAp = sin $ <5A M + cos $ <5Xm • (60) 

Unlike the modes discussed so far, the amplitude mode has finite eigenvalue in the x — limit. However, for x « 0.1 
the u> m = eigenvalue for the <j> mode has reached the value 0.493, quite close to the amplitude mode eigenvalue 
of 0.669. As seen from Eqs. (56,57,60) both these modes involve the fluctuation of the amplitudes of \ an d A, and 
they will admix for finite frequency. This complicates the interpretation of these collective modes. In contrast, the 
eigenvalue of the 9 mode is quite low at 0.155 and its interpretation as a collective mode is more clear. In Fig. 10 we 
show the spectral function for the ampltidue mode at q = (n, w). 




FIG. 10. Spectral function of the amplitude mode at q = (n, n). 

VI. EXPERIMENTAL OBSERVATION OF THE COLLECTIVE MODES 

Finally, we discuss possible experimental consequences of these collective modes. As discussed earlier, the 9 mode 
produces staggered orbital currents. These in turn generate a physical magnetic field which is staggered and which is, 
in principle, observable by inelastic neutron scattering. Here we compute the scattering cross-section. The neutron 
scattering by orbital currents has been considered by Hsu et al. 32 and we follow their discussion. We find that 



da k 



dfldq k 
where 

2 

i ] 

1 - e-P*" ] V* 



J (9' w ) = i \\ — \ (1- cos q x )lmxx,x(q,u) 



2 

q x 



+ y— J (1 - cos q y ) I mx y ,y(q,uj) 
-1(1 -e^) (l-e-^)lm Xx ,y(q,u;) 

- \ (1 - e-^) (1 - e***) ImXyAl,") ) (62) 
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where x^v — (J^Jv), A* = ~^-^^rzS ls the neutron magnetic moment, m n is the neutron mass, and V is the volume. 
We take the current operator to be the mean field expression (b = bo = y/x) 

J fl =ixt^2(jJ a f i+lJta -c.c^ . (63) 

icr 

In order to compute the correlation function, Xij.vi we note that any bilinear fermion operator can be written in 
terms of tpjAUijtpi for a suitable AU which in turn can be expanded according to Eq.(26) as AU^i +il — X^o 7 ?^ 7 " - 
The correlation function is then computed by treating rf^ as source terms and then differentiate the effective free 
energy with respect to 77°. The source terms simply modify the effective action S e g [Eq.(52)] by 

S' eS = X a M^Xp + X' a M^X> p (64) 

where X' a is obtained from X a by 6U" — ► 5U® + 77". Upon completing the square and integrating out X a , we find 
that (O^O* ) where 0° = i4+fj, T °" t f'i 1S obtained from the appropriate matrix element of 

M F (l-(M J + M F )- l M F ) . (65) 

In this way the neutron scattering cross-section is evaluated numericaly. We expect that the scattering is predomi- 
nantly coupled to the 9 mode and indeed the result is very similar to the 9 mode spectral function shown in Fig. 6. In 
Fig. 11 we show the lineshape I(q, to) given by Eq.(62) in the zero temperature limit at q — (tt, tt). In order to estimate 
the experimental feasibilty, we compare the total cross-section (integrated overall q and lo) to the scattering from a 
lattice of S = \ moments. Hsu et al. estimated the total cross-section to be about 1% of that of spin scattering and 
our results are in rough agreement. The experimental detection of this signal may be difficult because a resonance in 
the spin scattering exists around 30 meV at q — (tt, tt) for underdoped cuprates. This resonance is very narrow in q 
and u>, but its integrated weight near (tt, tt) is also about 1% of the total spin scattering. The scattering due to the 9 
mode is of comparable total strength but more spread out in lo, making it harder to detect. However, unlike the spin 
fluctuation which is isotropic, the orbital currents give rise to an effective moment which is perpendicular to the a-b 
plane. Since neutron is sensitive only to the comonent of the moment normal to the q vector, the orbital contribution 
can in principle be extracted by varying the q vector from normal to parallel to the a-b plane. 




FIG. 11. Neutron scattering intensity I(q,uj) at zero temperature vs. w (in units of J = |J) at q = (it, it). 

Another probe which couples to the orbital current is Raman scattering. However, since the excitation is epected 
to peak at (tt, tt), the large momentum transfer requires X-ray Raman scattering. It is known that X-ray couples to 
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magnetic moments as well as orbital moments. 33 The moment due to the orbital current is weak and corresponds 
in strength to roughly 0.1/!^. While Bragg scattering from an ordered spin moment has been observed, 34 inelastic 
scattering from a short-range order of such a small moment is beyond the current capability of X-ray scattering. We 
note, however, that polarization dependence is a powerful tood to distinguish between spin and orbital contributions 
in X-ray scattering. 33 

Next we discuss the possible measurement of the <j) gauge modes. As discussed before, these modes involve the 
modulation of the amplitude of the hopping parameter Xij an d the pairing parameter Ajj in a staggered fashion, i.e., 
with momentum (n, ir). The amplitude fluctuations of the superconducting order parameter is not easy to detect and 
has been experimentally observed only in the special case of superconductivity in a charge density wave system. 35:36 
On the other hand, the fluctuations in Xij couples to the quasiparticle hopping matrix element. In the mean field 
theory we can write an effective coupling as 

Hi = -xt ]T eW^flfj, + c.c. (66) 

<ij>0 

where we include Ais +Ix = aA^ as the lattice version of the electromagnetic field A. Expansion of Eq.(66) to first order 
in Aij yields the standard j ■ A coupling. When expanded to second order we obtain \xt (|) 2 A 2 - (^fj a fj a + c - c -) • 
This gives rise to Raman scattering which is coupled to fluctuations in Xij ■ This kind of coupling was discussed by 
Shastry and Shraiman 37 as an explanation of the continuum background due to incoherent electronic excitations. 
Here we expect that Raman scattering will couple to the transverse and longitudinal <fi gauge mode. Physically, a 
modulation of \xij\ is a modulation of the bond charge density which should couple to Raman scattering. Standard 
Raman scattering provides essentially zero momentum transfer. In order to couple to the 4> gauge mode at (ir, it) and 
to follow its dispersion, X-ray Raman scattering will be needed. The leading contribution to inelastic X-ray scattering 
originates from the 2 ^ ci A 2 term in the single particle Hamiltonian and the scattering cross-section is usually written 
as 38 



d 2 a ( do 



d^lduo \ dQ 



S(q,co) (67) 

Th 



where the Thomson cross-section is 



do~ \ 2 ( ~ » ,2 w 2 



- 



tq = e 2 /mc 2 is the Thomson radius, ei, u>\, and t2, ^2 arc the incident and scattering polarization vector and frequency, 
respectively, and S(q, u>) is the Fourier transform of the electron density-density correlation function. On a lattice, 
the corresponding matrix element comes from the second-order expansion of Eq.(66) and can be written as 

<i,j> 

-j(£)(s?)5?^<">(A'*~ + ") (69) 



where the sum is over nearest neighbors and ta — -r^r- Apart from which is of order ^ we see that the 



j 2 = t^-t- Apart from which is of order ^ 

4m* 1 m* 1 2 ' 

coupling is similar to the continuum theory, except that the number operator j ia ji a is replaced by \ ( f^fi+^a + ex.). 
Furthermore, a factor x arises due to the strong correlation. We see from Eq.(69) that X-ray Raman scattering directly 
couples to the fluctuating in Xij an d therefore to the cf) gauge mode. It seems that high resolution inelastic X-ray 
scattering is a promising technique to observe the appearance of the <j) gauge mode at low temperatures. The cross- 
section for X-ray scattering is proportional to < \-){q, uj)\ 2 > where 



20 



i(q, w) = ^2 ^x»(Q> ' ( 70 ) 

e and e' are the incoming and outgoing photon polarization and x M (q,tj) is the Fourier transform of the kinetic energy 
operator defined in Eq.(58). Thus the Raman scattering is given in terms of 

s* v (q>u)=<Xi*(<i>v)xv(q,u)> ■ (7i) 

This is computed numerically using the same method described by Eqs.(64) and (65). The results are shown in Fig. 
12 along (q x ,n). 




FIG. 12. Spectral function 5^(qr, w)(/U = x) for the correlation of the kinetic energy operator [Eq.(71)] which is measured 
by inelastic X-ray scattering. A contour plot and a bird's eye view are shown, q ranges from (0,tt) to (tv,tt), i.e., q — (q x ,n). 
Frequency uj is in units of J(= | J). 

The spectral function is peaked at (7r,7r) as expected. Note that along the (q^,7r) direction (Fig. 12) the spectral 
function has similar structure as the longitudinal </> gauge mode (Fig. 8). 

The spectral function at q = (tt,tt) is shown in Fig. 13. We note that its shape is quite similar to that of the 
amplitude mode shown in Fig. 10. For completeness, we also computed the correlation 5a =< SA^SA^ > where 
A M was defined in Eq.(59). This describes the fluctuations of the real part of Ajj and would correspond to the 
conventional superconducting amplitude mode. As shown in Fig. 14, its shape is similar to that of the <f> gauge mode. 
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FIG. 13. The lineshape of the kinetic operator correlator S* tli (q,u) (f/, = x) as a function of u) (in units of J = |J) at 
q = (7T, 7r). This is the curve labeled as RPA. Also shown is the curve labeled as bare bubble, which is computed by coupling 
to particle- hole excitations without allowing any collective enhancement. 




FIG. 14. The spectral function Sa(q,w) at q = (7r,7r). 



Given that S\ and SA are linear superpositions of SA and S<p according to Eqs. (56,60), it is surprising that their 
spectral functions are not simply weighted averages of the amplitude and (f> gauge modes. The reason is that there is 
significant cross-correlation between Sx and SA, as well as between SA and Sip, i.e., the amplitude and <fi fluctuations 
are not truly eigenmodes. This is a consequence of the significant mixing at x = 0.1 between these modes at finite 
frequencies as discussed earlier. It is perhaps better to focus on S% v which is experimentally measurable. The peak 
observed in corresponds to collective oscillations of the hopping matrix element. This is a new degree of freedom 
not present in conventional superconductors. To emphasize this point we have also computed the "bare bubble" 
response, i.e., using M£g without the collective enhancement shown in in Eq.(65). As seen in Fig. 13, the latter is 
much more spread out in frequency. Finally we note that for the relatively high doping considered here, the <f> gauge 
mode appears to be closely related to the more conventional superconductor amplitude fluctuation SA^ (albeit at 
q = (ir, 7r)) and is difficult to detect. 

We note that since Xij couples to the quasiparticle hopping matrix element, it in turn is coupled to phonons. In 
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Appendix B we consider the possibility that the collective modes may appear as phonon sidebands which may be 
detected optically. 

VII. CONCLUSION 

We have shown that in the gauge theory description of the t-J model, new collective modes appear in the su- 
perconducting ground state. The SU(2) gauge theory allows us to classify these collective modes and predict them 
analytically. These predictions are confirmed by numerical computation of the collective mode spectra. In particular, 
we see that the SU(2) formulation allows us to smoothly connect to the x = limit, whereas in the U(l) formulation, 
spurious modes would have appeared. 

In this paper we describe the bosons in the radial gauge. It is then clear that boson dynamics do not play a 
role and that collective modes dynamics are entirely determined by the fermions via the fermion bubbles which 
make up the matrix M F . In Ref. (24) we proposed a cr-model formulation where the low- lying excitations are 

parametrized by the boson h = Jxz = ^fx \ Zl ) . We proposed an effective action which included a term xz^D z 

where Do = d T + eAo + OqT 3 [Eq.(78) Ref. (24)]. Note that this term produces a coupling between the 8 and <f> 
modes of the form 88d T <fi. This Berry's phase term is absent in the present paper. Thus Eq.(78) of ref. (24) disagrees 
with the present work and we believe that it is incorrect. Within the cr-model approach in ref. (24), one should have 
absorbed the z'd T z term into the z^a^T^z term by a gauge transformation and solve for the a' terms by imposing 
the constraint locally. That is in fact what is done in this paper for small 89 and 8(f). A similar procedure can be 
adopted for arbitrary 8 and <j> to generate a cr-model. Thus the cr-model should be viewed as a way to parametrize the 
low- lying fluctuation of the Uij matrix [as in Eq.(40)] and the effective action depends only on the fermion dynamics. 

The appearance of the new collective modes answers the question of whether the superconducting state described 
by the gauge theory is any different from a conventional BCS state. The 8 mode is coupled to staggered orbital 
currents. The importance of these currents was already revealed in the Gutzwiller projected BCS wavefunction. 17 In 
principle, they should be observable by neutron scattering or X-ray Raman scattering. A new collective degree of 
freedom is the moduation of the hopping amplitude which is observable by high resolution inelastic X-ray scattering. 
We also show how they give rise to side-bands in certain phonon modes. However, the weakness of the coupling makes 
its observation difficult. The predictions of these collective modes are unique features of the slave-boson/gauge- field 
approach to the t-J model and experimental searches for the collective modes will serve as important tests of this line 
of approach to the high T c problem. 
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APPENDIX A: DERIVATION OF THE EFFECTIVE ACTION FOR COLLECTIVE MODES 

In this appendix, we sketch the calculation of the effective action for the collective modes up to quadratic order. 
For simplicity, we consider here the case of t' = 0. The inclusion of i'-tcrm is rather trivial. We start with the 
Lagrangian Eq.(21) in the text, because we consider the superconducting ground state with the bose condensation. 
The procedure is standard as in the usual 1 /N expansion. Namely we first divide the integral variables into the mean 
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field (saddle point) values and the fluctuations as given in Eq.(24). Integrating over the fermionic integral variables, 
we can expand the effective action with respect to the fluctuating part SUij, Sciq, 6oq, Sciq, and SR. 

Then the stationary condition that the linear order terms in 6Uij etc. vanish gives the self-consistent mean field 
equations as 

= [ n dkx T dk y -g fc 7 fc = r d 2 k -t^Tk ( 
X ~ ]_„ 2tt J_„ 2tt 2E k J (2tt) 2 2E k [ > 

J (2tt)2 E k 1 > 

Ao = i(2^^r (A3) 



where j k = cosk x + cosk y , (3 k = cosk x — cosk y , and E k = \J1$. + A'f, with £ k and A k being given below eq.(23). 
The Gaussian fluctuations are represented by the quadratic terms S2. 



>■ = ! dr 





J ]T (I^+mI 2 + l*A ii+ „| 2 ) - ti^2(2iS<4(t)6Ri + X (SRi) 2 ) 



■ ,H=x,y 



- r 2 t ij SR i 5R j <flf ja + f}j ia > 

<i,j>a 

+ ^Tr[GVGV] (A4) 

where f3 is the inverse temperature and the last term represents the second order contributions from the fermionic 
determinant. Here the Green's function is 2 x 2 matrix for each k = (k,u n ) (k: wave vector, Lo n : fermionic Matsubara 
frequency). Then the explicit form of the last term in eq.(A4) is (q = (q,u m )) 

2^ S E tr ^ G{k - i^ v ( k - ^ k + i/ 2 w k + ii^ vi - k + «/ 2 ' k - «/ 2 )]> ( A5 ) 

q k 

where tr means the trace over the 2x2 matrix. The Green's function is given by, 

which is represented in a compact form in terms of the Pauli matrices r 1 , r 2 , t 3 as 

G(k) = -K 2 + S 2 ]" 1 ^ + AfeT 1 - &T 3 ]. (A7) 

The matrix V is explicitely given by 

9 T 

V(k + q/2,k-q/2) = -= 

v f H=x,y 



-iSal(q) - Sx^(q)sink^ - 5 X '^{q) cos fc p , -i5a\{q) - Sal(q) + SA'^q) cos fc M + iSA'^q) cos k p 

-i5al(q) +Sal(q) + <5A^(q) cos fc p - iSA'^q) cos fc M , +^((7) - S X '^(q) sinfc M + ^(?) cos k^ 

Here it is noted that SUij is represented as 

5U ii+ll = -S x '^)t 3 + iSx'tt) + SA'^t 1 - SA'^t 2 (A9) 

where Xl ^(i) = XuW +»Xu(')i an d A M (i) = A' (i) +iA''(i). The Fourier transformations of these variables are defined 
as 



(A8) 
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(A10) 



and similar expressions for other variables. We define the 12 variables X a (q) (a = 1, 12) as 

1 5 X ';(q) x 

6A'(q) 



X(q) 



5A''(g) 
<%(?) 
Sx' v (q) 
SR(q) 
Sal(q) 
Sal(q) 
V S4(q) ) 

In the calculation of the fermionic polarization function the following integrals are needed. 

F ab (q, to m ; *0 = ^ 2 ^[^(fc + q/2)r a G(k - q/2)r b ], 



(All) 



(A12) 



where r° = 1 (unit matrix) and r 1 ' 2,3 are Pauli matrices. At zero temperature, the summation over Lo n is reduced to 
the integral, i.e., — > J" cLo/2w, which can be done to result in the following 4x4 matrix. 

F{q,w m ) = 

/ (w 2 j4 + x+X-+A + A_)7 1 -7 3 , i(A + + A_)7 2 +iu; m (A_ - A+)/i/2, 
z(A+ + A_)/ 2 + iw ro (A_ - A+)/ 1 /2, (u4/4 - x+X- + A+A_)/i - 7 3 , 

i( X+ A_ - X-A+)/i, (-x+ + X-)7> + ^ m (x+ + X-KiA 
V -»(X+ + X-)h + iuj m (x+ - X-)h/2, -(x+A_ + X-A+)7 1; 



(x+ - X-)h - 
( w m/ 4 - X+X 



where 

h{q) 



\ fill /VT 

(A+-A_)7 2 



^ m + 2col(El+E 2 _) + (E% -Elf 



A_ + X-A+)/i, -i(x+ + X-)h + ioJ m {x+ - X-)h/ 
2 - w m (x+ + X-)h/2, -(X+A_ + X-A+)/i 

_X_ - A+A_)7i - 7 3 , (-A+ + A_)7 2 + cj m (A+ + A_)Ji/ 
- w m (A + + A_)/i/2, + x+X- - A + A_)7i - 7 3 



(A13) 



col - El + El , o£ + 



2£ 4 



2E- 



h{q) ^ m + 2^ m {E\+El) + {E\-E 

1 



2 ^2 



UO 



hiq) < + 2col(El + El) + (El El 
The static limit of these funcitons are easily estimated as 



(£ + -£_)K 2 n + (£+-£-) 2 ] 

4E + E_ 

(£+ + £_)< | (E + +E^ujI | (E + -E_Y(E + +E_) 
SE-^-E— SE-^-E— 2 



(A14) 



^4^ = 0)= 2e+e _ {e++e _ } 

h{q,u m = o) = o 



(A15) 



Here we have 

the variables X a is given by 



introduced the abbrebiations such as E± = E k ± q / 2 etc. Then the quadratic action S 2 with respect to 
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Here 



q L Q ,/3=1,12 



^2 JX a {q)X a {-q) - ix(X 12 (q)X 9 (-q) + X 12 (q)X 9 {-q)) - ^-(cosq x + cosq y )X 9 (q)X 9 (-q). 



a=l,8 



n Q /3 = ^^2C a (k)Ci3(k)F a{a)Mfj) (q,uj m ;k) 



(A16) 



(A17) 



where 



—2 sin k x 
—2 sin k y 
2 cos k x 
2 cos k y 
—2 cos k x 
—2 cos k y 
—2 cos k x 
—2 cos fc„ 



V 



-Axi{cos(k x ) cos(q x /2) - 
-1 
-1 
-1 



cos(fc y ) cos(g y /2)) 



(A18) 



/ 



and a(a) is the index of the Pauli matrix corresponding to each component and is given as a(l) = a(2) = 0,a(3) = 
a(4) = 1, a(5) = a(6) = 2, a(7) = a(8) = a(9) = 3, a(10) = 1, a(ll) = 2, a(12) = 3. The fcrmion bubble is analytically 
continued in the standard way as F(iuj m lo + irj). Equation (A16) is numerically evaluated by discrctizing the 1st 
Brillouin zone by 200 x 200 and treating i] as small but finite. The convergence with respect to the number of the 
lattice points has been checked. 

Equation (A15) is the quadatic forms for the 12 variables X a , and by integrating over the last 3 variables 
X w , Xn, Xi 2 (= Sa 9 , 8oq, 8clq), we obtain the effective action for 9 variables, which has positive definite engenval- 
ues for each q = (q,uj m ) when the mean field solution is stable. 



APPENDIX B: COUPLING OF THE <j> GAUGE MODE TO PHONONS 



Let us consider the mode where the planar oxygen moves in and out of the plane. We focus on the oxygen mode 
because the mass is light and the frequency relatively low, and both features tend to enhance the coupling. In LSCO 
and YBCO, the oxygen is buckled out of the plane. (In YBCO the displacement mq = 0.256 Angstroms.) Then the 
out-of-plane phonon mode has a linear coupling to the Cu-0 bond length and therefore to the effective hopping t and 
exchange J. This problem was considered by Normand et al. 39 who concluded that in YBCO 

Stij = XftSuij/a (Bl) 
SJij = XjJSuij/a (B2) 

where they estimate X t w ^Aj w 2.6 for a displacement du^ of the oxygen on the ij bond normal to the plane. The 
surprising large coupling is partly due to the fact that the displacement Siiij is normalized to the lattice constant 
a w 4 Angstroms which is quite large. Let us first focus on the Stij term. The modulation of the Hamiltonian is 

SH t = ]T 8tijc\„Cj„ 

<ij>cr 
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<ij>,a 

nil ■ ■ 

*Xt(Xo + Sxij)xt-^ ( B3 ) 

The last line is in the mean field approximation, where we retained the amplitude fluctuation Sxij of \ij ■ 
Next, recall that the <j> gauge mode couples to \ij via Eqs.(37) and (38), so that 

S Xi j = -M-l) ix+iv Vij ■ (B4) 

Combining Eqs.(B3) and (B4), we find an effective coupling between the pnonon displacement and the <j> guage mode 
co-ordinate 

nil ■ ■ 

H eff = -X eff (-iy^^ Vij (B5) 

where 

Xeff = XtxtA Q . (B6) 

Note that a fluctuation of t>y couples to a (n, n) phonon mode as expected. It turns out that the modulation of J 
given by Eq.(B2) does not couple to the <f> gauge mode. The reason is that in mean field, the modulation of the J 
term is 8Jij(XijXij + Ay Ay). This couples to the amplitude mode. On the other hand, in the <\> gauge mode, the 
total modulus \Xij\ 2 + |Ay| 2 is held fixed, as shown in Fig. 2. Thus there is no coupling between the <j) gauge mode 
and the phonon via the (5 Jy term. 

We next approximate the phonon as an Einstein model. The energy is given by \(-u 2 + uj 2 )a 2 M(^f) where 
uj is the frequency and M is the mass. We approximate the gauge mode as a well defined mode at frequency iv^ 
with energy \^{—u 2 + w|)uy and adjust to match the spectral weight. The coupled phonon-0 gauge modes are 
obtained by diagonalizing the 2x2 matrix 

The modes are 



D -l = 2<V " '^1 (B7) 

V -Xeff \{-ui 2 +u 2 )a 2 M) ' 



1 



(co 2 + u %) ± ^( w g - u 2 ) 2 + Ag 2 ^ j (B8) 
where the dimensionless coupling constant is 



9 „ 2 U^a,, 

The phonon Green's function D = (SuijSuij) is 



{-uj 2 +uj 2 ) 

D = {-oj 2 +uj 2 + )(-LU 2 +LU 2 _)a 2 M ' (B10) 

Assuming uo > uj^, the phonon mode is shifted slightly upwards in frequency to uo + but a side band appears at cj_ 
which is near uj,. The ratio of the spectral weight of the w_ side band to the oj + mode is 

/ ,2 , ,2 

R=^ =-^± . (Bll) 

For weak coupling g 2 « 1, we expand the numerator using Eq.(B8) and replace w+ by oj n and cj_ by uj^ to obtain 
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This result for the relative spectral weight is also valid if ujq < uo<p. 

The spectral weight ratio is mainly determined by the coupling constant g 2 . To make a rough estimate, we take 
M = oxygen mass uj q w J/4 so that a 2 Muj 2 w 25 eV w 200 J. Taking A e // ~ J/4, we find 

2^J_^_ (B13) 

y 200w o> 

We estimate a^, w 1 because we can see from Fig. 9 that the spectral weight of the <j> gauge mode at (n, ir), is of order 
unity. Unfortunately, g 2 turns out to be very small. 

The small spectral weight (~ 5x 10~ 3 of the main phonon peak) means that it is probably impossible to observe the 
side-band by neutron scattering. Optical measurements offer higher precision. For YBCO the phonon is at q — (tt, n) 
and do not couple to light. LSCO offers a special opportunity in that at low temperatures, the lattice is in the low 
temperature orthorhombic phase (LTO) where the buckeling of the oxygen is staggered, i.e., Uij = uq{— 1)*=+% . Then 
a q = distortion Uij = + Suij couples to the staggered modulation in Hence in LSCO the (j> gauge mode 
is coupled to the q = phonon mode which can be studied optically either by absorption or by Raman, depending 
on its activity. The signature of the side band is that it should appear only in the superconducting phase, because 
in the normal state (pseudogap state) the I vector rotates out of the plane and is disordered, so that the <j> mode is 
expected to be smeared out. Similarly, the weight of the side band will be reduced by applying a magnetic field by 
the fraction of the sample occupied by the vortex core. This is because the I vector is rotated to the north pole near 
the vortex core and the <f> mode will lose its identity. 

The above discussion was based on the assumption that the <j) gauge mode is a well defined mode and that the 
amplitude mode SA is much higher in frequency and plays no role. This is correct for small doping but we have 
seen that for x = 0.1, the numerical results show a strong admixture of the <f> gauge mode and the amplitude mode. 
Nevertheless, it is still correct that the phonon couples to , which produces a sideband with a lineshape given 
in Fig. 13. However we should include the fact that the phonon also modulates the exchange constant according to 
Eq.(B2) and couples to the amplitude fluctuation as well. Since the amplitude fluctuation has a very similar lineshape 
(see Fig. 10) to that of S^, the final result may still be interpreted as a modulation of the hopping amplitude, a new 
collective degree of freedom not present in conventional superconductors. Finally, we remark that we only have results 
for T = 0, and in the x = 0.1 case the issue of how much of the spectral weight in S*^ survives above T c remains 
open. 
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